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TECHNICAL PUBLICATION 


AN EXAMINATION OF SELECTED GEOMAGNETIC INDICES IN RELATION 

TO THE SUNSPOT CYCLE 

1. INTRODUCTION 


Over the years, a number of geomagnetic indices have been devised to describe the variation of 
the geomagnetic field in response to the changing conditions of the solar wind at Earth. 1 - 2 A few of these 
indices, in fact, have been shown to reliably predict the strength of the following sunspot cycle several 
years in advance. 3-14 It is this aspect that will be examined more closely in this Technical Publication. 

This study is divided into seven results sections. Section 2.1 examines the observed yearly average 
values of the aa and Ap geomagnetic indices, as well as the number of disturbed days (NDD), defined as 
those days when the Ap index equals or exceeds 25, all in relation to annual sunspot number ( R ). The aa 
index was introduced by Mayaud in the early 1970s 15 - 16 and is defined as the average of 3-hr K indices 
converted to the amplitude of the field at two nearly antipodal observatories (in England and Australia). 
The aa index, as derived by Mayaud, is available from 1868 to the present. Although Nevanlinna and 
Kataja 17 have provided an extension for the aa index back to 1844 using Helsinki magnetic records, these 
extended values have not been used in the present study. Ap, or the daily equivalent planetary amplitude, 
is a daily index of geomagnetic activity that is determined using a linear scale rather than the quasi- 
logarithmic scale of the K indices. Values for it and NDD are available from 1932 to the present. (Yearly 
values of R, aa, and Ap can be found online at ftp://ftp.ngdc.noaa. gov/STP/ . 18 ) 

Section 2.2 examines the effects of suggested changes to the observed aa record. Svalgaard, 
Cliver, and Le Sager 19 derived a new geomagnetic index called the Inter- Hour Variability (IHV) index 
for investigations of the long-term variability of the solar wind-magnetosphere system and used it to 
reconstruct the observed aa record. For the interval 1957-2000, they found yearly averages of their proxy 
aa index (based on observations at Cheltenham/Fredericksburg during the two 3-hr periods between 
0 and 6 hr UT) and the IHV index to be virtually the same, while for the interval preceding 1957, there was 
considerable difference between the indices. Hence, they concluded that values of the observed aa index 
prior to 1957 might be inaccurate. 

Section 2.3 examines annual minimum values of aa, Ap, and NDD (aamm, Apmin, and NDDmin, 
respectively) for each cycle in relation to //max (the annual maximum amplitude of the sunspot cycle), 
both as single-variate and bivariate fits (the latter fit also using //min, the annual minimum amplitude of 
the sunspot cycle). Single-variate and bivariate fits previously have been shown to reliably predict the size 
of the ongoing sunspot cycle some 2 to 3 years in advance. 5 


Section 2.4 examines cyclic averages (minimum-to-minimum based on R ) of aa, Ap, NDD, and 
NSSC (the number of sudden storm commencements), denoted <aa>, <Ap>, <NDD>, and <NSSC>, 
respectively, and the cyclic sums (the total number recorded over the sunspot cycle) of NDD and NSSC, 
denoted 2NDD and 2NSSC, respectively, all in relation to the cyclic average of R, denoted <R>, and 
i?min and i?max. 

Section 2.5 examines Apmax and the number of days when Ap equals or exceeds 100 and 200, 
denoted ND(Ap> 100) and ND(Ap>200), respectively. 

Section 2.6 identifies single-variate and bivariate regressions of Ap, NDD, and Aprnax against 
aa and aa and R, respectively, for the interval of 1932-2005. These regressions are then used to provide 
estimates of Ap, NDD, and Apmax prior to 1932. 

Section 2.7 looks at epoch analyses of sunspot number and the various geomagnetic parameters 
(aa, Ap, NDD, A/; max, and NSSC) relative to E(Rmdx) and £(i?min) to determine when onset for cycle 24 
should be expected. 
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2. RESULTS AND DISCUSSION 


2.1 Geomagnetic Indices: aa,Ap, and NDD 

Figure 1 displays the annual variation of sunspot number R in (panel (a)), aa (panel (b)), Ap 
(panel (c)), and NDD (panel (d)). For R and aa, values are given for 1868-2005, while for Ap and NDD, 
values are given for 1932-2005. The lower numbers refer to sunspot cycles 11-23. The thin vertical lines 
refer to the occurrences of the sunspot minimum year for each of the sunspot cycles. 

Inspection of figure 1 reveals two important aspects of solar/geomagnetic activity: First, there has 
been a general rise in R and aa over the years and second, while R usually is single peaked, aa (and Ap and 
NDD) generally has two or more peaks, with the first being closely associated with the rising/maximum 
phase of the sunspot cycle and the other(s) occurring during the declining portion of the sunspot cycle, 
often occurring just prior to the onset of the new cycle. Hathaway, Wilson, and Reichmann 20 previously 
have shown that the secular trend in R is long-term, extending from the Maunder minimum. Also, Clilverd 
et al., 21 using two long-running European stations (Sodankyla and Niemegk), have reconstructed the long- 
term aa index and concluded that the trend in aa is real, inferring a long-term increase in solar coronal 
magnetic field strength. Indeed, Wilson and Hathaway 22 have reported that, on the basis of 10-yr moving 
averages, R and aa are highly correlated (r=0.933). 

In order to explain the later-occurring peaks in the aa index (which nearly always have been the 
largest of the cycle), Feynman 23 suggested that two components comprise the aa index — one that is in 
phase and correlated directly with the sunspot cycle and the other (the residual) that is out of phase and 
associated with interplanetary disturbances from the Sun (high-speed solar wind streams and the like 24 ). 
Figure 2 plots aa versus R (panel (a)), Ap versus R (panel (b)), and NDD versus R (panel (c)). The identified 
lines in each panel are those that allow for a determination of the inferred sunspot cycle component for 
the geomagnetic indices by passing a line through the two lowest points. Hence, all geomagnetic index 
values lie either on or above these lines. By subtracting this component from the observed value of the 
geomagnetic index, one infers the value of the residual interplanetary component. (The technique employed 
here is slightly different from that employed in Hathaway, Wilson, and Reichmann 13 and Hathaway and 
Wilson. 14 ) 
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Figure 1. Variation of annual averages of sunspot number, R (panel (a)); aa-geomagnetic index, 
aa (panel (b)); Ap-geomagnetic index, Ap (panel (c)); and number of disturbed days, 
NDD (panel (d)). See text for details. 
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(c) 0 100 200 


Annual Mean Sunspot Number, R 

Figure 2. Scatterplots of aa versus R (panel (a)); Ap versus R (panel (b)); 
and NDD versus R (panel (c)). All values of the geomagnetic 
indices lie on or above the identified regression lines. 
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Figure 3 plots the residual interplanetary components of acij (panel (a)), Ap l (panel (b)), and NDDy 
(panel (c)), all plotted as the thin jagged lines. The thicker, smoother lines are 2-yr moving averages 25 of 
the residuals (using a weighting of 1 :2: 1), and, as before, the lower numbers refer to the individual sunspot 
cycles 11-23 and the thin vertical lines mark the occurrences of the sunspot minimum years. 





Figure 3. Variation of the interplanetary components aa l (panel (a)); Ap l (panel (b)); 
and NDD / (panel (c)). The thick lines are 2-yr moving averages of the 
residual interplanetary components. 
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Figure 4 shows histograms marking the frequency of occurrences for the maximum amplitude 
of aci/ (E(aa, max, panel (a)) and for the last occurrences of the local maximums of acij (E(aa, max (last), 
panel (b)), which often also turns out to be E{aa l max) for most of the cycles (11, 12, 14, 15, 17, 18, 21, and 
22), in years relative to the occurrences of EfEmax). For E{aa I max), two cycles (1 1 and 16) have E(aa /max) 
2 yr after EfEmax); four cycles (12, 13, 19, and 23) have E(aajmax) 3 yr after E(Rmax); five cycles (14, 15, 
18, 21, and 22) have E(aajmax) 5 yr after £(Emax); and two cycles (17 and 20) have E(aajmax) 6 yr after 
E^Emax). For E(aajmax (last)), the bulk of the cycles (13, 14, 15, 18, 21, 22, and 23) have E(aajmax (last)) 
5 yr after E(Rmax). 
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Figure 4. Histograms of the frequency of occurrences for Eiaapnax) 
(panel (a)) and E((aapnax (last)) (panel (b)), relative to the 
elapsed time in years from E(Rmax). 


Figure 5 displays scatterplots of Emin (cycle n+ 1) versus acij max (cycle n) (panel (a)); Emin 
(cycle /t+l) versus Apjmax (cycle n ) (panel (b)); Emin (cycle n + 1) versus NDD 7 max (cycle n) 
(panel (c)); Emax (cycle n+ 1) versus aa : max (cycle n) (panel (d)); Emax (cycle n+ 1) versus Ap l max 
(cycle n) (panel (e)); and Emax (cycle n + 1) versus NDDpnax (cycle n) (panel (f)). In each panel, the 
small, downward-pointing arrow marks the parametric value for cycle 23. For panels (b), (c), (e), and (f), 
no statistically significant (confidence level ( cl) >95 percent) linear correlation can be inferred between 
the parameters. For these, only the mean and standard deviation ( sd ) are shown. However, for panels (a) 
and (d), statistically significant positive linear correlations are inferred (those related to aaj max). Thus, 
given aa / max = 26.1 for cycle 23, one infers Emin= 13 ±4.9 and Emax = 183.7 ±46.3 for cycle 24, the next 
sunspot cycle. Both predictions are 90-percent prediction intervals so, there is only a 5-percent chance that 
Emin will be smaller than 9.1 and Emax will be smaller than 137.4. Likewise, there is only a 5-percent 
chance that Emin will be larger than 17.9 and Emax will be larger than 230. 
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Figure 5. Scatterplots of AAnin (cycle n+ 1) versus selected maximum interplanetary 
components for cycle n (panels (a)-(c)) and A 1 max (cycle n+ 1) versus 
selected maximum interplanetary components for cycle n (panels (d)-(f)). 
See text for detail. 


Figure 6 depicts scatterplots of 7?min (cycle n + 1) versus aa / max(last) (cycle n) (panel (a)); /?min 
(cycle n + 1) versus Ap 7 max(last) (cycle n) (panel (b)); i?min (cycle n + 1) versus NDD / max (last) (cycle n) 
(panel (c)); A 1 max (cycle n + 1) versus acijmax (last) (cycle n) (panel (d)); Ahnax (cycle n + 1) versus Ap 7 max 
(last) (cycle n) (panel (e)); and i?max (cycle n+ 1) versus NDD 7 max (last) (cycle n) (panel (f)). For panels 
(b) and (c), no statistically significant linear correlation can be inferred between the parameters. For these, 
only the mean and sd are shown. For the other panels, statistically significant positive linear correlations 
are inferred. Thus, given aa 7 max (last) = 15 for cycle 23 (an assumed value, since the values for 2006 and 
2007 are presently unknown, true also for Apjmax (last) and NDD/inax (last)), one infers Ahum = 7.2 ±5 
and /?max= 1 17.5 ±41.2 for cycle 24. Also, given Apjmax (last) = 5.8 and NDD 7 max (last) = 3 1 for cycle 
23, one infers A!max= 110.1 ±38.4 and A’max = 96.2 ±35.6, respectively, for cycle 24. These predictions 
clearly differ with the aforementioned predictions based on acijmax for cycle 23. 

Figure 7 shows scatterplots of 7?min 9 (cycle n + 1) versus act/ max 2 (cycle n) (panel (a)), where 
the subscript 2 refers to the 2-yr moving average (shown in fig. 3 as the thicker smoother lines); /?min 2 
(cycle n + 1) versus A/; / max 2 (cycle n) (panel (b)); A?min 2 (cycle n + 1) versus NDD / max 2 (panel (c)); 
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Figure 6. Scatterplots of /On in (cycle n+ 1) versus selected last-occurring 
maximum interplanetary components for cycle n (panels (a)-(c)) 
and //max (cycle n+\) versus selected last-occurring maximum 
interplanetary components for cycle n (panels (d)-(f)). 

See text for details. 


//max 2 (cycle n + 1) versus aa / max 2 (cycle n) (panel (d)); //max 2 (cycle n+\) versus Ap l max 2 (cycle n ) 
(panel (e)); and //max 2 (cycle n+ 1) versus NDD 7 max 2 (cycle n ) (panel (f)). Only panels (a) and (d) are 
found to contain a statistically significant positive linear correlation, both having higher coefficients of 
correlation (r) and smaller standard errors (se) than for either of the unsmoothed plots (figs. 5 and 6). Thus, 
given aa / max 2 =18.6 for cycle 23, one infers //min 2 = 16.9±4.8 and //max 2 = 146.7 ±30.5 for cycle 24, 
both predictions being 90-percent prediction intervals. Hence, there is only a 5-percent chance that cycle 
24’s /On in 2 will be smaller than 12.1 or larger than 21.7, and there is only a 5-percent chance that 
cycle 24’s //max 2 will be smaller than 116.2 or larger than 177.2. The 90-percent prediction interval for cycle 
24’s //min 2 suggests that its /On in (the annual average value) will be about 10.4 ±3 and its Rm (12-mo 
moving average, or smoothed monthly mean sunspot number minimum amplitude) will be about 9.4 ± 2.9. 
Similarly, the 90-percent prediction interval for cycle 24’s /Zmax 2 suggests that its //max (the annual 
average value) will be about 157.5 ±31.8 and its RM (12-mo moving average, or smoothed monthly mean 
sunspot number maximum amplitude) will be about 162.9 ±33.5. Thus, these values suggest that cycle 
24’s activity will be greater than average (and larger than was seen for cycle 23), similar in size to that 
experienced in cycles 21 and 22. 14 ’ 26 ’ 27 
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Figure 7. Scatterplots of /?min 2 (cycle n+ 1) versus selected 2-yr moving averages of the 
maximum interplanetary components for cycle n (panels (a)-(c)) and i?max 2 
(cycle 7t+l) versus 2-yr moving averages of the maximum interplanetary 
components for cycle n (panels (d)-(f)). See text for details. 


2.2 aa (adjusted) Geomagnetic Index 

Recall from figure 1 the long-term rise in the observed aa index over time, which is also apparent 
in R, this indicating a strong correlation between the parameters. Indeed, as previously mentioned, Wilson 
and Hathaway 2 - have noted that on the basis of 10-yr moving averages, aa and R are highly correlated 
(r=0.933). Also, as previously mentioned, Svalgaard, Cliver, and Le Sager 19 in their reconstruction of 
the observed aa index using the IHV index found that, for the interval of 1957-2000, both parameters are 
virtually identical, while for the interval preceding 1957, the reconstructed aa (based on the IHV index) 
lies above the observed aa. Hence, they concluded that the observed aa prior to 1957 might be inaccurate. 
(In 1957, the Northern Hemisphere (NH) magnetometer used in deriving the aa index was moved from 
Abinger, England, to its present location in Hartland, England.) 
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Figure 8. Variation of selected parametric differences, 1901-2000. 
See text for details. 


Figure 8 shows the differences aa— ^chl/FRD > (P an el (a)); < 3 < 3 chl/FRD _ ^V, A 2 (panel (b)); 
and aa- IHV, A 3 (panel (c)), where aa is the observed aa index value; ««chl/FRD ‘ s the derived aa value 
based on observations at Cheltenham/Fredericksburg during the two 3-hr periods between 0 and 6 hr UT 
(which represented the observed or proxy aa index in the Svalgaard, Cliver, and Le Sager study 19 ), and 
IHV is the reconstructed aa value, these latter two values taken from Svalgaard, Cliver, and Le Sager. 
Across the top of the chart are time ticks 1-5, corresponding to the notes given to the right of the figure. 
For example, note 1 says that values of IHV prior to 1915 are reduced 30 percent, as advised in Svalgaard, 
Cliver, and Le Sager; note 2 states that the Southern Hemisphere (SH) magnetometer was moved from 
Melbourne, Australia, to Toolangui, Australia, and so forth. Today, the NH magnetometer is located at 
Hartland, England, and the SH magnetometer is located at Canberra, Australia. These stations are now 
determining the official aa index. 

Inspection of figure 8 reveals that each of the differences appears to be relatively stable from 1957 
onward, averaging about 2.3 for Aj, 0. 1 for A ? , and 2.3 for A 3 . For the earlier portion of the record, 1901— 
1956, Aj is found to average about 1.8, A 2 about -3.1, and A 3 about -1.3. The differences of the means, 
especially for A ? and A 3 , are found to be statistically significant. Hence, an offset-correction appears to be 
needed for the earlier interval (1901-1956), one being about 3 for A 2 and about 3.6 for A 3 . 
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Figure 9 repeats the analysis, but this time comparing aa, Ap, and IHV. Plotted are the differences 
aa-Ap (A 4 ) and Ap-IHV (A g ). The notes at the top of the chart refer to the notes explained in figure 8. For 
the interval 1957-2005, A 4 averages about 9.1 and A 5 averages about -6.7, while for the earlier interval 
of 1932-1956, they average, respectively, 6.8 and -7.5. So, it appears that, relative to Ap, the aa index 
needs to be offset-corrected by 2.3 units and, relative to IHV, aa should be offset-corrected by about 3.1 
units (=2.3 + 0.8). 


4 5 




Figure 9. Variation of selected parametric differences, 1932-2000. 
See text for details. 


Figure 10 replots aa versus R, but this time adding in an offset value of 3 to the observed values 
of aa for all values of aa prior to 1957. In truth, while slight offset-corrections might be attributable after 
each movement of the NH or SH magnetometers, for simplicity sake, a uniform correction of 3 is used to 
account for the general difference inferred between the earlier (prior to 1957) and later aa datasets. 
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Figure 10. Scatterplot of aa (adjusted) versus R. 


Figure 11 replots the aaj component of the adjusted aa dataset, following the format used in figure 3. 
While the values of at// (adjusted) are shifted upward, as compared to that shown for aaj in figure 3, there 
remains noticeable in both figures an increase in the value of aa. For example, for the interval 1868-1956, 
the observed (unadjusted) aa averaged 17.17, while for the interval 1957-2005, it averaged 23.98, an 
increase of nearly 40 percent. Using the adjusted aa, for the first interval, it averaged 20.17, as compared 
to 23.98 for the second interval, an increase of about 19 percent. For the same two intervals, aa R averaged 
9.69 and 12.06, nearly a 25-percent increase; aa R (adjusted) averaged 11.79 and 13.61, an increase of about 
15 percent; aaj averaged 7.48 and 11.92, an increase of about 59 percent; and aaj (adjusted) averaged 8.38 
and 10.37, an increase of about 24 percent. Thus, the increase in geomagnetic activity appears real, being 
seen in both the sunspot cycle component and the interplanetary component, regardless of whether the 
data have been adjusted or not. 
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Figure 11. Variation of aa^adjusted). The thick line is the 2-yr moving 
average of the residual interplanetary component. 


Figure 12 plots i?min 2 (cycle n + 1) versus at// (adjusted) max 2 (cycle n) (panel (a)), where the latter 
parameter refers to the maximum 2-yr moving average of the aa l (adjusted) as plotted in figure 1 1 , for each 
cycle; and i?max 2 (cycle n + 1) versus aa / (adjusted)max 2 (cycle n) (panel (b)). The downward-pointing 
arrow gives the value of aa 7 (adjusted) max 2 for cycle 23, equal to 16.9. Thus, for cycle 24, the 90-percent 
prediction interval for 7?min 2 = 16±4.7 and i?max 2 = 143.2 ±23.2. These values are slightly lower than 
what was gleaned using the observed aa (unadjusted) values. In terms of i?min and Rm, the 90-percent 
prediction intervals for cycle 24 are 9.8 ±2.9 and 8.8 ±2.8. Similarly, in terms of i?max and RM , the 90- 
percent prediction intervals for cycle 24 are 153.8 ±24.7 and 159 ±25.5. 
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Figure 12. Scatterplots of i?min 2 (cycle n+ 1) versus aa / (adjusted)max 2 (cycle n ) (panel (a)) and 
tRmax 2 (cycle n + 1) versus acij( adj u s ted ) m a x 2 (cycle n) (panel (b)). See text for details. 
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2.3 Single- Variate and Bivariate Fits Based on Parametric Minimum Values 

Figure 13 displays cyclic values of //min (panel (a)); A 1 max (panel (b)); aamin (panel (c)); Apmin 
(panel (d)); and NDDmin (panel (e)) for cycles 12-23. For aamin, Apmin, and NDDmin, these are the 
minimum values found in the vicinity of sunspot minimum (within 1 yr following sunspot minimum). This 
caveat is necessary because cycle 21 actually had lower values in 1980 during the maximum amplitude 
phase of the sunspot cycle, the only cycle in the entire record to have such an anomaly. Also, for aamin, 
two lines are plotted: the values based on the observed record and the values based on the adjusted record. 
Additionally, statistically significant regressions are given for //min, /An ax, and aamin (both observed and 
adjusted). For Apmin and NDDmin, only the mean and sd are given. 

Figure 13 clearly shows that over the course of cycles 12-23, there has been an apparent secular rise 
in //min, /An ax, and aamin. On the basis of these inferred trends, one estimates the 90-percent prediction 
interval for each of the parameters to be //min= 12.3 ±4.8 (or Rm= 1 1 .2 ± 4.6), //max= 167.9±55 (or 
RM= 173.3 ±56.5), aamin (observed) = 20.6 ±4.6, and aamin (adjusted) = 20.4 ±5.2. 

Figure 14 displays single-variate scatterplots of /An ax versus /On in (panel (a)); aamin (panel (b)); 
A/;m in (panel (c)); and NDDmin (panel (d)). Thus, observation of the minimum values of these parameters 
allows (at least for aamin, Apmin, and NDDmin) for the prediction of //max for the ongoing sunspot cycle 
some 2-3 yr in advance. (Although neither plot of //max versus //min is statistically significant, that based 
on cycles 12-23 or that based on cycles 17-23, the regressions are shown because they will be used in the 
next figure as part of a bivariate fit.) 

From figure 14, it should be obvious that, if cycle 24 has an aamin of about 20, as suggested by 
the linear secular trend, then, clearly, one should expect an //max of about 160 or so for cycle 24. Such a 
value also suggests that Apmin should measure about 11 and NDDmin about 25 or so. Through July 2006, 
Ap has averaged only about 7, compared to 13.6 for the year 2005, and the total number of disturbed days 
amounts to only 9, compared to 43 for the year 2005. Hence, either these values will substantially increase 
as the year progresses, especially in the fall of the year, or cycle 24 ’s values will fall well below the 
means. It is interesting to note that based on a presumed repeating pattern of “below-above-above-below” 
observed in Apmin and NDDmin, one expects both Apmin and NDDmin to be “above” their respective 
means, suggesting that the latter half of 2006 might be geophysically quite active. 
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Figure 13. Cyclic variation of selected parameters. Notice the strong 
secular increases in i?min, /driax, and aa. 


17 


200 i 


yi 2—23 = 77.033 + 5.521 x 
r= 0.525, r 2 = 0.276 
se = 36.320, c/>90% 


200 “I 


y = 14.206 + 7.686X 
r= 0.890, r 2 = 0.793 
se = 19.427, cl >99.9% 


y' (adjusted) 

: y 



(a) 0 


Rmin 


10 


15 


y' = -32.343 + 9.736x 
r= 0.940, r 2 = 0.883 
se = 14.588, c/ >99. 




>Apmin 



Figure 14. Scatterplots of jRmax versus selected minimum values of the parameters. 


Figure 15 shows scatterplots of i?max (observed) versus idnax (predicted) on the basis of specific 
bivariate fits, where bivariate fit 1 ( bv j) is based on /dnin and amnin (panel (a)); bv 2 on 7?min and 
aamin (adjusted) (panel (b)); bv 2 on /dnin and Apmin (panel (c)); and bv 4 on tRmin and NDDmin (panel 
(d)). For the various fits, the one having the smallest standard error of estimate is bv 4 . Thus, once /dnin and 
NDDmin are observed for cycle 24 (perhaps in 2007 or 2008), an alternate means for predicting Rm&x for 
cycle 24 will be available (having a 90-percent prediction interval measuring about ±15.5 units of sunspot 
number). 
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Figure 15. Scatterplots of Rnrdx (observed) versus selected bivariate fits. 
See text for details. 


2.4 Cyclic Averages and Sums 

Figure 16 displays cyclic averages and sums, from sunspot minimum to sunspot minimum, of 
<R> (panel (a)); <aa>, both the observed and adjusted (panel (b)), where the dashed line represents 
the adjusted values; <Ap> (panel (c)); <NDD> (panel (d)), where the dashed line refers to 2NDD; and 
<NSSC> (panel (e)), where the dashed line refers to 2NSSC. (<R> for cycle 11 excludes the year 1867 
value of R and all parameters for cycle 23 exclude the year 2006 values.) 

For <R>, <aa>, <aa> (adjusted), <NSSC>, and 2NSSC, all display statistically significant 
upward linear secular trends, given by the regression equations. Thus, on the basis of the inferred trends, 
the 90-percent prediction interval of the parameters for cycle 24 is <R> = 82.4 ±29.5, <aa> = 26.5±4, 
<aa> (adjusted) = 26.2 ±4.2, <NSSC>=35.2±6.1, and 2NSSC = 358±62. Also, because of the presumed 
“below-above-above -below” pattern, one expects <Ap> to be above 14.8, <NDD> to be above 54.4, and 
2NDD to be above 563. 
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Figure 16. Variation of cyclic averages of selected parameters. 
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Figures 17 and 18 show comparisons of each of the parameters against i?min (left panels) and 
against i?max (right panels). All statistically significant regressions (c/> 95 percent) are identified. Hence, 
once i?min has been observed for cycle 24, one can use these regressions to estimate the cyclic parametric 
averages (and sums). Likewise, once Rma.x has occurred, one can use these regressions to estimate with 
greater accuracy the cyclic parametric averages (and sums). (One can also use the secular estimates to 
estimate i?min and jRmax for cycle 24. The secular estimate for <R> suggests that i?min will be about 14 
and /x’max about 160; the secular estimate for <aa>, whether adjusted or not, suggests i?min will be about 
15 and i?max about 185; the secular estimate for <NSSC> suggests i?min will be about 15 and tRmax about 
190; and the secular estimate for 2NSSC suggests i?min will be about 15 and i?max about 200.) 

A peculiarity noted in figure 18 is the cyclic split seen in <NSSC> and 2NSSC. Namely, cycles 
11-16 and 17-23 appear to be grouped distinctly from each other. Such groupings reinforce the notion 
that cycles of late have been more robust as compared to earlier cycles. If the trend continues, then 
cycle 24, likewise, should be expected to be another robust cycle. 26 ’ 28-30 (Another peculiarity is the 
apparent “below-above-above-below” pattern noticeable in <Ap>, <NDD>, and 2NDD versus i?min or 
tRmax in figs. 17 and 18, which, if real, suggests above average values for cycle 24.) 

Figure 19 shows the scatterplot of yearly counts of NSSC versus R. Plainly, years of higher R are 
associated with years of higher NSSC, and vice versa, although even during periods of low R (near sunspot 
minimum), one still expects to see about 10 or more sudden storm commencements. 
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Figure 17. Scatterplots of cyclic averages of selected parameters against i?min (left panels) and 
7?max (right panels). The numbered filled circles refer to individual sunspot cycles. 
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Figure 18. Scatterplots of additional cyclic averages and sums 
against i?min (left panels) and A 1 max (right panels). 
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Figure 19. Scatterplot of NSSC versus R. 
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2.5 Apmax, ND ( Ap > 100) and ND (Ap > 200) 

Figure 20 displays the yearly Apmax for the interval 1932-2005. As before, the numbers at the 
bottom of the chart refer to sunspot cycles 17-23 and the thin vertical lines refer to the epochs of sunspot 
minimum for each cycle. Noticeable is that all cycles have two or more peaks of activity, with the largest 
Apmax values usually occurring after i?max (two exceptions: cycles 18 and 22). Also, all cycles have 
peak values in excess of Ap= 100 and all but one (cycle 20) had peak values exceeding Ap = 200, with the 
highest Apmax (=280) occurring in cycle 19 in November 1960 (after its i?max). For cycle 23, its highest 
Apmax occurred in October 2003 (= 204), the sixth highest value on record. 

Figure 21 depicts the scatterplot of Apmax versus R. Although there is a statistically significant 
positive correlation between the parameters, associating higher (lower) Apmax with higher (lower) R, one 
clearly sees that even when the sunspot cycle is near solar minimum there is still opportunity for large Ap. 
For example, for the 15 yr when R <20, Apmax spanned 38 (cycle 23) to 202 (cycle 22), with one-third of 
the years having Apmax> 100. 
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Figure 20. Variation of the annual maximum daily Ap value, Apmax. 



Figure 21. Scatterplot of Apmax versus R. 
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Figure 22 displays the scatterplot of Apmax versus Ap. Plainly, higher (lower) A/; max associates 
with higher (lower) Ap and the association is stronger (r = 0.67) than is found for R (r=0.54). For Ap < 10, 
as yet, there has never been an A/; max > 100. (82 has been the largest value of Apmax for Ap< 10.) 



Figure 22. Scatterplot of A/; max versus Ap. 


Figure 23 shows the frequency of occurrence for E(Apmax), defined here as the epoch of A/; max 
occurrence, during the year. Apparent is a semi-annual variation of the geomagnetic field, with two 
prominent peaks near the equinoxes and minima near the solstices. 31 

Figure 24 gives the frequency of occurrence of E(Apmax) relative to E(Rmax), defined here as the 
epoch of sunspot maximum occurrence, for cycles 17-23. Most of the cycles have their Apmax value 3 to 
4 years after Efftmax), as has previously been noted. 

Figure 25 displays the cyclic variation of A/; max (panel (a)), and ND (Ap > 100) and ND (Ap > 200) 
(panel (b)) for cycles 17-23. Clearly, cycle 19 has been the most geomagnetically active cycle on record, 
having the highest A/unax (=280) and the largest ND(A/;> 100) and ND(A/;>200) (39 days and 5 days, 
respectively). All cycles have had at least 13 days (cycle 20) with Ap > 100, with cycle 23 experiencing 22 
such days so far, including one day of Ap > 200. 
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Figure 23. Semi-annual variation of Apmax. 
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Figure 24. Histogram of the frequency of occurrences of E(Apmax) 
relative to the elapsed time in years relative to £(7?max). 
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Figure 25. Cyclic variation of Apmax (panel (a)) and ND(A/;> 100) 
and ND ( Ap > 200) (panel (b)). 
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Figure 26 depicts scatterplots of Apmax versus ft min (panel (a)); A/; max versus ftmax (panel (b)); 
ND (Ap> 100) versus ftmin (panel (c)); and ND(Ap> 100) versus ftmax (panel (d)). Of the plots, only 
the latter is found to show a statistically significant positive linear regression. Thus, given the size of an 
ongoing sunspot cycle, one can crudely estimate the number of days during that cycle when Ap> 100. 
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Figure 26. Scatterplots of A/; max versus ftmin and ftmax (panels (a) and (b), 
respectively) and ND(Ap> 100) versus ftmin and ftmax 
(panels (c) and (d), respectively). 


2.6 Estimating Ap, NDD, and Apmax Prior to 1932 

Tables 1 and 2 identify single- variate (against aa or aa (adjusted)) and bivariate (ft and act or 
aa (adjusted)) regressions, respectively, for Ap, NDD, and Apmax for the interval 1932-2005. All 
regressions are highly statistically significant and might prove useful for determining estimates of Ap, 
NDD, and Apmax for the earlier interval 1868-1931. It should be noted that the bivariate fits offer no 
significant improvement over the single-variate fits based on aa or aa (adjusted) alone. 
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Table 1. Selected single-variate regressions for 1932-2005 (ft = 74). 


Correlation 

r 


se 

c/(%) 

Ap = -2.906 + 0.761 aa 

0.947 

0.896 

1.32 

>99.9 

Ap = -4.806 + 0.808aa (adjusted) 

0.970 

0.942 

0.99 

>99.9 

NDD = -63.375 + 5.078aa 

0.926 

0.857 

10.32 

>99.9 

NDD = -77.194 + 5.438aa (adjusted) 

0.957 

0.916 

7.91 

>99.9 

Apmax = -37.966 + 7. 1 21 aa 

0.657 

0.431 

40.59 

>99.9 

Apmax = -47.403 + 7.21 3aa (adjusted) 

0.642 

0.412 

41.26 

>99.9 


Table 2. Selected bivariate regressions for 1932-2005 (ft = 74). 


Correlation 

R Y.12 

R2 Y.12 

S Y.12 

Ap = -2.684 + 0.733aa + 0.006R 

0.948 

0.900 

1.27 

Ap = -4.51 0 + 0.775aa (adjusted) + 0.007R 

0.973 

0.996 

0.94 

NDD = -64.855 + 5.279aa-0.044R 

0.929 

0.863 

10.15 

NDD = -78.434 + 5.600aa (adjusted) - 0.037R 

0.959 

0.920 

7.73 

Apmax = -27.488 + 5. 636aa + 0.332R 

0.716 

0.512 

37.85 

Apmax = -35.376 + 5.667 aa (adjusted) + 0.352R 

0.711 

0.506 

38.11 


Table 3 provides the estimates of Ap, NDD, and Apmax for 1868-1931 using the aforementioned 
regression fits. For simplicity, the tabular entries for NDD and Apmax have been rounded to the nearest 
whole number. Each parametric column contains two numbers: The first is the value based on the single- 
variate fit and the second is based on the bivariate fit. As an example, for 1868, the entries for Ap are 10.9 
(single- variate fit) and 10.9 (bivariate fit); for Ap (adjusted), they are 12.3 (single- variate fit) and 12.2 
(bivariate fit), and so on. To obtain the 90-percent prediction interval for the estimate, it is approximately 
1.668 times the appropriate se for the single- variate fits and 1.668 times the appropriate S Y i 2 f° r the 
bivariate fits (given in tables 1 and 2). Thus, for Ap, the 1868 estimate is 10.9 ±2.2 (single- variate fit) 
and 10. 9± 1.5 (bivariate fit); for Ap(adjusted) it is 12.3 ± 1.7 (single-variate fit) and 12.2 ± 1.6, and so on. 
(The negative numbers in table 3 for NDD should be interpreted as 0+ 1.668 se or 0+ 1.668 5^12’ where 
se equals 10.3 or 7.9 for the single-variate fits, dependent upon using the observed or adjusted data, 
respectively; S Y 12 equals 10.2 or 7.7 for the bivariate fits, dependent upon using the observed or adjusted 
data. 
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Table 3. Parametric estimates for 1868-1931. 


Year 

Ap 

Ap ( adjusted ) 

NDD 

NDD ( adjusted ) 

Apmax 

Apmax ( adjusted ) 

1868 

10.9 

10.9 

12.3 

12.2 

29 

30 

38 

39 

92 

88 

106 

81 

1869 

12.9 

13.0 

14.4 

14.5 

42 

42 

52 

52 

110 

114 

124 

126 

1870 

14.0 

14.4 

15.6 

16.0 

49 

46 

60 

57 

120 

144 

134 

156 

1871 

13.3 

13.6 

14.8 

15.1 

45 

43 

55 

53 

114 

130 

128 

142 

1872 

15.1 

15.2 

16.7 

16.8 

57 

55 

68 

67 

130 

139 

145 

151 

1873 

12.5 

12.5 

13.9 

13.9 

39 

39 

49 

49 

106 

108 

120 

119 

1874 

8.2 

8.3 

9.4 

9.4 

11 

10 

19 

18 

67 

70 

80 

80 

1875 

5.6 

5.6 

6.7 

6.6 

-7 

-7 

0 

0 

42 

41 

55 

51 

1876 

4.3 

4.3 

5.3 

5.3 

-15 

-15 

-9 

-9 

30 

30 

43 

39 

1877 

3.9 

3.9 

4.8 

4.8 

-18 

-18 

-13 

-13 

25 

27 

38 

36 

1878 

2.6 

2.6 

3.4 

3.4 

-27 

-27 

-22 

-22 

13 

14 

26 

24 

1879 

2.4 

2.5 

3.3 

3.3 

-28 

-28 

-23 

-23 

12 

14 

25 

23 

1880 

5.8 

5.9 

6.9 

7.0 

-5 

-6 

2 

1 

44 

48 

57 

58 

1881 

7.4 

7.6 

8.6 

8.7 

6 

5 

13 

12 

59 

67 

80 

78 

1882 

14.5 

14.5 

16.1 

16.0 

53 

53 

64 

64 

125 

121 

139 

132 

1883 

10.4 

10.5 

11.8 

11.8 

26 

25 

34 

34 

87 

92 

101 

103 

1884 

7.8 

8.0 

9.0 

9.2 

8 

7 

16 

15 

62 

73 

76 

84 

1885 

8.8 

8.9 

10.1 

10.1 

15 

14 

23 

22 

72 

77 

85 

87 

1886 

12.8 

12.6 

14.3 

14.0 

41 

43 

51 

53 

109 

97 

123 

107 

1887 

9.6 

9.4 

10.9 

10.6 

20 

21 

28 

29 

79 

69 

93 

79 

1888 

8.8 

8.6 

10.1 

9.8 

15 

16 

23 

24 

72 

62 

85 

71 

1889 

6.6 

6.5 

7.7 

7.5 

0 

1 

7 

8 

51 

45 

64 

55 

1890 

5.2 

5.2 

6.3 

6.2 

-9 

-9 

-3 

-2 

38 

35 

51 

45 

1891 

10.0 

10.0 

11.4 

11.2 

23 

23 

32 

32 

83 

80 

97 

91 

1892 

15.5 

15.5 

17.2 

17.1 

60 

60 

71 

71 

134 

133 

149 

145 

1893 

10.0 

10.3 

11.4 

11.6 

23 

21 

32 

30 

83 

97 

97 

108 

1894 

12.8 

13.0 

14.3 

14.4 

42 

41 

52 

51 

109 

115 

124 

126 

1895 

10.9 

11.0 

12.2 

12.3 

29 

28 

38 

37 

91 

96 

105 

107 

1896 

10.7 

10.7 

12.1 

12.0 

28 

28 

37 

37 

90 

88 

103 

98 

1897 

7.4 

7.4 

8.5 

8.5 

5 

5 

13 

13 

58 

57 

72 

67 

1898 

8.6 

8.5 

9.8 

9.7 

13 

14 

21 

22 

70 

67 

83 

77 

1899 

7.1 

7.0 

8.2 

8.1 

3 

3 

10 

11 

55 

50 

69 

60 

1900 

2.8 

2.9 

3.7 

3.7 

-25 

-26 

-20 

-20 

15 

18 

28 

28 

1901 

1.7 

1.7 

2.5 

2.5 

-33 

-33 

-28 

-28 

5 

7 

18 

17 

1902 

2.0 

2.1 

2.9 

2.9 

-30 

-31 

-26 

-26 

8 

11 

22 

20 

1903 

6.1 

6.2 

7.2 

7.2 

-3 

-3 

4 

4 

47 

48 

60 

58 

1904 

5.9 

6.1 

7.0 

7.1 

-5 

-6 

2 

2 

45 

52 

58 

62 

1905 

8.4 

8.6 

9.7 

9.8 

12 

11 

20 

19 

68 

78 

82 

88 

1906 

6.5 

6.7 

7.6 

7.8 

0 

-2 

7 

6 

50 

60 

64 

71 

1907 

9.3 

9.4 

10.5 

10.6 

18 

17 

26 

25 

76 

83 

90 

94 

1908 

10.0 

10.1 

11.4 

11.3 

23 

23 

32 

32 

83 

84 

97 

95 

1909 

10.1 

10.1 

11.4 

11.4 

24 

24 

32 

32 

84 

84 

98 

94 
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Table 3. Parametric estimates for 1868-1931 (Continued). 


Year 

Ap 

Ap (adjusted) 

NDD 

NDD (adjusted) 

Apmax 

Apmax (adjusted) 

1910 

10.4 

10.3 

11.8 

11.5 

26 

27 

34 

35 

87 

77 

101 

87 

1911 

9.1 

8.9 

10.4 

10.1 

17 

18 

25 

26 

75 

64 

88 

73 

1912 

3.8 

3.8 

4.7 

4.7 

-19 

-19 

-13 

-13 

25 

23 

38 

33 

1913 

3.6 

3.6 

4.6 

4.5 

-20 

-20 

-14 

-14 

23 

21 

36 

31 

1914 

5.4 

5.4 

6.4 

6.3 

-8 

-8 

-2 

-1 

40 

37 

53 

47 

1915 

9.0 

9.0 

10.2 

10.2 

16 

15 

24 

24 

73 

76 

87 

87 

1916 

12.2 

12.2 

13.6 

13.6 

37 

37 

47 

47 

103 

103 

117 

114 

1917 

10.9 

11.3 

12.3 

12.6 

29 

27 

38 

36 

92 

110 

106 

121 

1918 

13.5 

13.6 

15.0 

15.0 

46 

45 

56 

56 

115 

120 

129 

132 

1919 

14.1 

14.1 

15.7 

15.6 

50 

51 

61 

61 

122 

120 

136 

131 

1920 

10.4 

10.4 

11.8 

11.6 

26 

26 

34 

35 

87 

84 

101 

94 

1921 

9.7 

9.6 

11.0 

10.8 

20 

21 

29 

30 

80 

74 

93 

84 

1922 

11.3 

11.1 

12.7 

12.4 

32 

33 

41 

42 

95 

83 

109 

93 

1923 

4.9 

4.8 

5.9 

5.8 

-12 

-11 

-5 

-5 

35 

32 

48 

42 

1924 

4.8 

4.8 

5.8 

5.8 

-12 

-12 

-6 

-6 

34 

35 

47 

45 

1925 

7.0 

7.1 

8.1 

8.2 

3 

2 

10 

9 

55 

61 

68 

71 

1926 

12.2 

12.2 

13.6 

13.6 

37 

37 

47 

47 

103 

105 

117 

116 

1927 

9.7 

9.9 

11.0 

11.2 

21 

20 

29 

29 

80 

89 

94 

100 

1928 

10.5 

10.7 

11.8 

12.0 

26 

25 

35 

34 

87 

98 

101 

82 

1929 

11.8 

11.9 

13.2 

13.2 

35 

34 

44 

44 

100 

103 

113 

114 

1930 

18.8 

18.4 

20.6 

20.2 

81 

84 

94 

96 

165 

145 

180 

156 

1931 

9.8 

9.7 

11.1 

10.9 

21 

22 

30 

31 

81 

74 

95 

84 


2.7 Epoch Analyses 

Figure 27 shows the comparisons between yearly averages of various parametric values for cycle 23 
(the filled circles) and the mean of the parametric values (the line) for cycles 11-22 (for R, aa, and NSSC) 
or 16-22 ( Ap , NDD, and Apmax) for elapsed time in years from E(R max) on the basis of epoch analyses. 
For R (panel (a)), cycle 23 appears to closely mimic the behavior of the mean of cycles 11-22, suggesting, 
perhaps, that minimum amplitude for cycle 24 (its onset) will occur in year 7 following cycle 23 ’s E(Rmax), 
this year corresponding to 2007. The proportion of cycles having £’(i?max) rt + j sooner than year 7 is 3/12, 
while the proportion having E{Rmax) n + \ in year 7 is 5/12 and the proportion having E(Rmax) n + i in year 
8 is 4/12. All cycles having E(Rmax) n + j in year 8 are cycles of longer period (minimum-to-minimum 
duration >135 mo on the basis of 12-mo moving averages or smoothed monthly mean sunspot number), 
while all cycles having E(Rmax) n+ j in year 7 are cycles of shorter period (minimum-to-minimum duration 
<126 mo). Those having E(Rmax) n + j in years 5 or 6 are a mixed bag, with cycles 15 and 16 being cycles 
of shorter period and cycle 12 being of longer period. Now, midway through 2006 (year 6), R has averaged 
only about 16. This value should become smaller as the year progresses, although it already is smaller 
than was seen in 5 of the past 6 cycles for comparison year 6 following E(R max). Also, it is nearly within 
the 90-percent prediction interval for //min based on cycles 11-22 (i?min 90 = 7±7). So, presently, one 
cannot discount that year 2006 might ultimately be the minimum amplitude year for cycle 24 marking its 
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Figure 27. Comparison of cycle 23 parametric values (filled circles) and cyclic averages 
(cycles 11-22 or 16-22) based on epoch analyses using //(//max) as the 
temporal marker. See text for details. 


onset. (It should be noted that, as yet, no new cycle spots have been seen, this being a crucial parameter 
for heralding the onset of a new sunspot cycle. 32-34 ) 


For aa (panel (b)), it too has mimicked the mean for cycles 11-22 rather closely, except for the 
value at year 3 from E(Rmax) which is the highest on record, being nearly 2.5 standard deviations above 
the mean. The bulk of the cycles (6/12) have had E(aamm ) n+ ( in year 8, corresponding to 2008, suggesting 
that minimum amplitude occurrence for cycle 24 should be expected in 2007, since E{aaram) usually 
follows F(//min) by 1 yr (true for 9 of 12 cycles). 
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For NSSC (panel (c)), again, it too mimics the mean for cycles 11-22, except now the value for 
year 3 from E(Rmax) is lower than any of the past 6 cycles, but within the observed range of values 
(9-49) for all known cycles. Flalf of the sunspot cycles had E(NSSCmin) /!+ ] in year 6 from EfRmax), 
corresponding to 2006. Since NSSC is strongly correlated (r =0.87) against R (fig. 19), it may be that 2006 
might turn out to be the minimum amplitude year for cycle 24. (It should be noted that cycle 21 had three 
consecutive years of NSSC =20, in years 5, 6, and 7 from E(Rmax). Here, year 6 has been used as the 
minimum value year.) 

For Ap (panel (d)), cycle 23 ’s values generally have been below that of the mean for cycles 16-22, 
except for the peak in year 3 from E(Rmax). However, whereas the peak at year 3 for cut is of record value, 
the value at year 3 for Ap is only the third largest, below that which was seen for cycles 19 (23.7) and 21 
( 22 . 6 ). 


For NDD (panel (e)), like Ap, the values generally have been below the mean for cycle 16-22, 
except for year 3 from Fornax). The value of 113 disturbed days, however, matches the highest ever seen, 
having occurred twice before in cycle 18 during years 4 and 5 from £’(i?max). 

For Apmax (panel (f)), cycle 23 ’s values have been above the mean for every year except year 2 
from E(Rmax). Its maximum value (Apmax = 204) occurred in year 3. 

Figure 28 is similar to figure 27, except now the comparison epoch is £’(7?min) w + 1 , meaning the 
epoch of the next cycle’s minimum amplitude 7?min. The filled circles refer to cycle 23 ’s parametric values, 
drawn presuming that onset for cycle 24 will be 2007, and the lines represent the parametric means. For R 
(panel (a)), the 90-percent prediction interval for R in year -1 relative to £(i?min) w+ j is /? 90 = 13.9 ± 11.1 
and, as already been mentioned, for year 0 relative to E(Rmin) n+l it is 7 ± 7. (If eventually the year 2006 is 
recognized as the minimum amplitude year for cycle 24, then cycle 23 ’s parametric values must be moved 
1 yr to the right in all panels.) 

For aa (panel (b)), yearly values for cycle 23 are expected to continue to decline through E(Rm'm) n + , 
(actually to the year following the onset year for cycle 24). Certainly, if 2007 represents the minimum 
amplitude year for cycle 24, then E(aamm) n+ j for cycle 24 should be expected in 2008. (The proportion 
of cycles having E{aamm) in the year following £’(7?min) is 9/12.) 

For NSSC (panel (c)), all yearly values for cycle 23 are above the mean, except for year-4 (2003). 
Its value (NSSC= 11) is more than two standard deviations below the mean and is outside the range of 
previously observed cycles, 14-49. On the other hand, one should note the sharp decrease in the mean 
at year -3. Should the low value of NSSC for the year 2003 be better associated with year -3 rather than 
-4, this might be an indication that the minimum amplitude year for cycle 24 is 2006 rather than 2007. 
Assuming that the 2003 value of NSSC is year -3, its value is now only one standard deviation below 
the mean and is within the range of previously observed cycles, 9-37. The proportion of cycles having 
£’(NSSCmin) w+ j at year-1 is 4/12 and 5/12 at year 0. (If one assumes the 2003 value represents the 
year -3 value from E(Rmm) n+ 1, then the last filled circle, the 2005 value, occurs at year-1 from 
£(7?min) w+1 .) 
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For Ap (panel (d)), all yearly values for cycle 23 are below the mean, except for year-4 (2003). 
Its value (Ap = 21.7) is within the range of previously observed cycles, whether one presumes it to be 
year -4 (12.7-23.7) or year -3 (14.4-22.3). The proportion of cycles having £’(Apmin) w+1 at year 1 
relative to E(Rmin) n+l is 5/7. 

For NDD (panel (e)), all yearly values for cycle 23 are below the mean, except for year-4 (2003), 
which is outside the range of previously observed cycles (32-107). Its value (NDD= 113), however, is 
within the range of previously observed cycles for years -3 (52-113) and -2 (40-113). The proportion of 
cycles having F(NDDmin) n+ , at year 1 relative to ii(i?min) n+ j is 7/7. 

For Apmax (panel (f)), all yearly values for cycle 23 have been above the mean, except for year -5 
(2002), although it is within the range of previously observed cycles (73-236). Shifting the 2002 value (78) 
to year-4, however, places it outside the range of previously observed cycles (136-280). The proportion 
of cycles having £’((Apmax)min) w+ j is 3/7 for year 0 and 4/7 for year 1 relative to E(Rm\n) n+ \. 






Elapsed Time in Years From E(Rmin)„ + 1 Elapsed Time in Years From E(/?min)„ + 1 Elapsed Time in Years From E(Rmin)„ + 1 


Figure 28. Comparison of cycle 23 parametric values (filled circles) 
and cyclic averages (cycles 11-22 or 16-22) based on 
epoch analyses using E(Rmax) n+l . See text for details. 
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3. CONCLUSION 


This study has examined a variety of geomagnetic indices in relation to the sizes and occurrences 
of minimum and maximum amplitudes of sunspot cycles. The study confirms that both sunspot number 
and the aa geomagnetic index have increased over time, such that cycles of late are among the largest on 
record. While both sunspot number and the aa index is strongly correlated, there are differences between 
them. Namely, each sunspot cycle usually has a single peak, while the aa index usually has two or more 
peaks, often with the major peak occurring during the declining portion of the sunspot cycle. 

The aa index was decomposed into two components — one associated directly with sunspot number 
R and the other being the residual interplanetary component due to coronal holes and the like. Plots of 
the residual interplanetary component, however, continue to display multiple peaks during the declining 
portion of the sunspot cycle, with the largest usually being the last peak just prior to sunspot minimum for 
the next sunspot cycle. The peaks, whether using the largest or the last-occurring peak, strongly correlate 
with both the minimum and maximum amplitudes of the following cycle, typically several years in 
advance. Unfortunately, the multiple peaks during the decline of cycle 23 present a dilemma for estimating 
the maximum amplitude of the next cycle. Based on the maximum value of the residual for cycle 23, 
cycle 24 will have 7?min 90 = 13 ±4.9 and /?max 90 = 183.7 ±46.3, while based on the last observed residual 
peak, cycle 24 will have i?min 90 = 7.2±5 and i?max 90 = 1 1 7.5 ±41 .2, yielding an overlap of 8.1-12.2 for 
tRmin and 137.4-158.7 for i?max. The dilemma can be mitigated by applying a 2-yr moving average to the 
residuals, thereby yielding predictions of i?min 90 = 10.4±3 and i?max 90 = 157.5 ±31.8 for cycle 24, these 
annual averages equivalent to 9.4±2.9 and 162.9±33.5, respectively, expressed as smoothed monthly 
mean sunspot number. 

In recent years, while studies maintain that the aa index has increased substantially over time, 
comparisons of it with the IHV index suggest that values of the aa index prior to 1957 might be somewhat 
inaccurate. In 1957, the NH magnetometer used for deriving the official aa index was moved from 
Abinger, England, to its present location in Hartland, England. The offset for the earlier aa values appears 
to measure about 3. Hence, by simply adding 3 to the observed values of aa, one can adjust aa, bringing 
the two disparate datasets into closer agreement. Doing so and repeating the analyses results in predictions 
of i?min 90 = 9.8±2.9 and i?max 90 = 153.8 ±24.7 for cycle 24, these values equivalent to 8.8±2.8 and 
159 ± 25.5, respectively, expressed as smoothed monthly mean sunspot number. Using the adjusted values 
for aa does not significantly alter the predictions of minimum and maximum amplitude for cycle 24, 
although it does reduce the uncertainty by 3 percent for i?min and 22 percent for /v’max. 

Once parametric cyclic minimum values are observed (usually, the minimum values of the 
geomagnetic indices occur in the year following sunspot minimum), one can use these minimum values 
of the geomagnetic indices to deduce with higher precision (±15.5 units of sunspot number) the size 
for the ongoing cycle, either on the basis of single-variate or bivariate fits (where the bivariate fits also 
incorporate i?min). If the epoch of sunspot minimum is the year 2006, one should expect the minimum 
in the geomagnetic indices in 2007. On the other hand, if the epoch of sunspot minimum is delayed until 
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2007, then the minimum in the geomagnetic indices should not be expected until 2008, with sunspot 
maximum amplitude expected about 2-3 yr later. The best predictor of sunspot maximum amplitude 
after the onset of the new cycle appears to be the one incorporating both NDDmin and ftmin (bv 4 ), being 
ftmax 90 = 121.4-5.877 ftmin + 3.665 NDDmin±15.5 and having r=0.978 and r 2 =0.956, meaning that 
95.6 percent of the variance can be explained by the bivariate fit. 

Each of the parameters ftmin, ftmax, and aamin displays a secular rise over cycles 12-23, 
such that values for cycle 24 are ftmin 90 = 12.3 ±4.8 (or ftm 90 = 1 1 .2 ± 4.6), ftmax 90 = 167.9 ±55 (or 
RMg 0 = 173.3 ±56.5), and aa 90 = 20.6 ±4.6. Secular rises in Ap, NDD, or Apmax cannot be detected, owing 
to the brevity of these records (only about half as long as that for aa ). 

Figure 29 shows the values for seven cycle-related parameters for the interval of January 2005- 
July 2006, including R, aa, Ap, NDD, Apmax, NSD (the number of spotless days 28,30 ’ 35 during the 
month), and Jmax (the maximum daily value of sunspot number during the month) (panels (a) through 
(g), respectively). Plotted are the monthly means of each parameter, with the horizontal line drawn during 
the year 2005 representing the yearly parametric average. To the right are yearly parametric averages 
(90-percent probability intervals) during the sunspot minimum year. Clearly, values now being experienced 
in the year 2006 suggest that cycle minimum is imminent; the yearly averages for 2006 are now falling 
within the 90-percent probability intervals indicative of the sunspot minimum year. 

In addition to the secular rises found for ftmin, ftmax, and aamin, cyclic averages (minimum-to- 
minimum) for some of the parameters are found to show statistically significant secular increases as well. 
For cycle 24, <ft> 90 = 82.4 ±29.5, <aa> 90 = 26.5 ±4, <NSSC> 90 = 35.2 ±6.1, and 2NSSC 90 = 35 8 ±62. 
Also, based on a presumed inherent behavior (“below-above-above-below”), <Ap> will be >14.8, <NDD> 
> 54.4, and 2NDD >563 for cycle 24. Furthermore, cyclic averages for cycle 24 can be deduced from 
cycle 24’s ftmin and ftmax, once they are observed. An unexpected finding is that <NSSC> and 2NSSC 
scatterplots cluster into two distinct cycle groupings: 11-16 and 17-23. 

Concerning A/; max, all cycles are found to have had at least 1 day of A/;max> 100 and all but 
one (cycle 20) had at least 1 day of A/;max>200. Opportunities for large Apmax remain throughout the 
solar cycle, including near sunspot minimum. For the 15 yr when R was <20, A/; max spanned 38 to 202, 
with one-third of the years having at least 1 day of Apmax> 100. A/nnax > 100 has never been seen when 
Ap< 10. Once ft max for a sunspot cycle has been observed, the number of days that A/nnax > 1 00 can be 
easily estimated. 

On the basis of the statistically significant regressions between Ap, NDD, and Apmax against aa 
during the interval 1932-2005, estimates can be made for the earlier interval 1868-1931. 

Epoch analyses for each of the parameters suggest that sunspot minimum year for cycle 24 probably 
is the year 2007, although one cannot, as yet, rule out the year 2006 as the sunspot minimum year. For 
some of the parameters, a better fit is found using the year 2006 as the sunspot minimum year, while for 
others, the better fit is found using the year 2007 as the sunspot minimum year. 


37 



7 ± 6 . 6 * 


14.5 ± 9 . 9 * 


11 ± 2 . 7 * 


2.5 ± 1 . 2 * 



39.7 ± 22 . 2 * 


17 ± 10 * 



Figure 29. Parametric values for 2005 and 2006. 
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